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PACS 03.75.Ss - Degenerate Fermi gases 

Abstract -We study the problem of a mobile impurity of mass M interacting via a s-wave broad 
or narrow Feshbach resonance with a Fermi sea of particles of mass m. Truncating the Hilbert 
space to at most one pair of particle-hole excitations of the Fermi sea, we determine ground state 
properties of the polaronic branch other than its energy, namely the polaron quasiparticle residue 
Z, and the impurity-to-fermion pair correlation function G(x). We show that G(x) deviates from 
unity at large distances as — (A4 + B4 cos 2kpx) /(kyx) 4 , where kp is the Fermi momentum; since 
A4 > and B4 > 0, the polaron has a diverging rms radius and exhibits Friedel-like oscillations. In 
the weakly attractive limit, we obtain analytical results, that in particular detect the failure of the 
Hilbert space truncation for a diverging mass impurity, as expected from Anderson orthogonality 
catastrophe; at distances between ~ 1/fe and the asymptotic distance where the 1/x 4 law applies, 
they reveal that G(x) exhibits an intriguing multiscale structure. 



Introduction. — The physics of atomic Fermi gases 
'"Q lias recently experienced a fast development, thanks to the 
Ch Feshbach resonance technique that allows to tune the s- 



wave scattering length a of the interaction, and to obtain 
ighly degenerate strongly interacting Fermi gases [TJ[2]- 
The first realizations of spin polarized configurations [3J 
' '[|] asked for theoretical interpretation of the experimental 
results, and shortly afterwards it was proposed that, in 
£■ — .the strongly polarized case, the minority atoms dressed 
t-H by the Fermi sea of the majority atoms form a normal gas 
OO of quasiparticles called polarons [5j|6] , which agrees with 
f^*) the experimental phase diagram [31[7] . 
,— ' The basic single-polaron properties, such as its binding 
energy to the Fermi sea and its effective mass, are now well 
• .understood, both for broad [5,6,8,9 and narrow [10 12 
. £h f esnDacn resonances. Here we study the polaron ground 
state properties going beyond the energy search, for both 
5_i types of resonance. We study the quasiparticle residue Z, 
.^already investigated experimentally in P21Q3] and theo- 
retically in [TTJ[Tn], and the intriguing issue of the density 
distribution that surrounds the impurity, as characterized 
by the pair correlation function: This gives access to the 
spatial extension of the polaron, of potential important 
consequences on the properties of the polaronic gas. 

Model. — We consider in three dimensions an ideal 
gas of TV same-spin-state fermions of mass m, enclosed in a 
cubic quantization volume V with periodic boundary con- 
ditions. The gas is perturbed by an impurity, that is a dis- 



tinguishable particle, of mass M. The impurity interacts 
with each fermion resonantly on a s-wave broad or narrow 
Feshbach resonance, as described by a two-channel model 
[TbTUSj . with a Hamiltonian H written in [TU]: The parti- 
cles exist either in the form of fermions or impurity in the 
open channel, or in the form of a tightly bound fermion- 
to-impurity molecule in the closed channel. These two 
forms are coherently interconverted by the interchannel 
coupling of amplitude A 0. We restrict to the zero-range 
(or infinite momentum cut-off) limit, so that the inter- 
action is characterized by the s-wave scattering length a 
and the non- negative Feshbach length i?» [TO] . In terms of 
the effective coupling constant g and of the interchannel 
coupling A, one has 



2irh 2 



and i?* = 



nh 4 
AV' 



(1) 



where fi = mM / (m + M) is the reduced mass. A broad 
Feshbach resonance corresponds to i?* = 0. 

Polaronic ansatz. — Whereas the ground state of 
the system presents two branches, a polaronic branch and 
a dimeronic branch [8l rT0l - fT2lfT5ll20H21] , we restrict to the 
polaronic branch. The ground-state of the fermions 
is the usual Fermi sea ("FS"), of energy E-ps(N). We 
determine the ground state of a single impurity interacting 
with the N fermions under the approximation proposed in 



x We can neglect the open-channel fermionic interaction if its 
(background) scattering length a^g is M the van der Waals length. 
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[3] (and generalized to encompass the narrow resonance 
case in PHHH]), that truncates the Hilbert space to at 
most one pair of particle-hole excitations of the Fermi sea. 
For a zero total momentum, this corresponds to the ansatz 



M/m = 0.1505 



M/m= 1 



M/m = 6.6439 



k.q 



IPS), 
(2) 



where u 



' and 6 k are the creation operators of an im- 



purity, a fermion and a closed-channel molecule of wave 
vector k. The prime above the summation symbol means 
that the sum is restricted to q belonging to the Fermi sea 
of N fermions, and to k not belonging to that Fermi sea. 
The successive terms in (0) correspond in that order to 
the ones generated by repeated action of the Hamiltonian 
H on the A = polaronic ground state. One then has to 
minimize the expectation value of H within the ansatz @, 
with respect to the variational parameters 4>, <^ q and </>kq, 
with the constraint (V'poilV'poi) = 1. Expressing 0kq in 
' q in terms of <j>, as in [TU], one is left with 



terms of (/> q and 



a scalar implicit equation for the polaron energy counted 
with respect to Efs(N); in the thermodynamic limit: 



Ai?poi = Ep i — Efs(N) 



' d 3 q 1 

(2^Pq 



(3) 



where the prime on the integral over q means that it is 
restricted to the Fermi sea q < kp, with the Fermi mo- 
mentum &f related as usual to the mean density of the 
Fermi sea p = N/V by fcp = (6tt 2 p) 1 ^ 3 ■ The function of 
the energy in the denominator of the integrand is 



= 1 _ pk F p 2 R* 
q g n 2 h 2 TrH* V 
" d 3 k' ( 1 



(A£; po i 



2p 

(271)3 V^q-k' + E V - £q - AE pol ~ Wk^ 



(4) 



where ew = for the fermions, Eu = for the 

impurity, and the prime on the integral over k' means 
that it is restricted to k' > fcp. 



Quasiparticle residue. — The polaron is a well- 
defined quasiparticle if it has a non-zero quasiparticle 
residue Z, which is defined in the Green's function for- 
malism from the long imaginary-time decay of the Green's 
function [5]. Within the polaronic ansatz @, it was shown 
in [15] that simply Z = \4>\ 2 - We first write the amplitude 
</> q in terms of the denominator (j4]) as [10] : 



A_ 



(5) 



where A is a normalization factor. Then using Eq. ([5]) and 
the coupled equations for <p, q and </>kq (TO] , we get the 
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Fig. 1: (Color online) Quasiparticle residue ^ for various 
mass ratios: M/m = 0.1505 (column a), M/m — 1 (column 
b), M/m = 6.6439 (column c), and various Feshbach lengths: 
fc F i?« = (row 1), k F R, = 1 (row 2), k F R* = 10 (row 3). 
Black line (solid for Z > 1/2, dashed for Z < 1/2): numerical 
solution. Red line: Second order weakly attractive expansion 
(1221 1 , that diverges for R, > when s of Eq. (|19p tends to unity, 
and for R, = when l/(fcpa) — ¥ 0. Points in (c2): Experimen- 
tal data of [12 for k F R t ~ 0.9474. Vertical dotted lines [ID] : 
Polaron-to-dimeron crossing point [on the left (resp. right) of 
this line, the ground state is polaronic (resp. dimeronic)]. 



following expression for the residue: 

d 3 q 



<" d 3 kd 3 q 



1 

A 2 



1 



(2tt)3 V\ 



(27r) 6 V^q-k+ek-Eq- AE po 



2 -1 



(6) 



In Fig.[T]we plot Z as a function of 1/fcpa for various mass 
ratios M/m and reduced Feshbach lengths kpR*. We find 
that Z tends to 1 when a — > - , as expected, and to 
when a — > + . The polaronic ansatz a priori makes 
sense when Z is close to unity, and its accuracy becomes 
questionable when Z — > 0. In Fig. [I] we thus have plotted 
Z in dashed line for a predicted value below 1/2. In the 
weakly attractive limit, we shall give below a systematic 
expansion of Z up to second order in fcpa. 

Pair correlation function. — The pair correlation 
function G(x„ — x^) is proportional to the probability den- 
sity of finding a fermion at position x u knowing that the 
impurity is localized at x^. In terms of the fermionic and 
impurity field operators ^(xu) and ipd{*-d)'- 



G(x u - x d ) 



( 4>t ( x « ) 4>d ) ipd (xd)V>„ (x u ) ) 



PPd 



(7) 



Here p = N/V is the unperturbed mean fermionic den- 
sity, and pd = (ip d (xd)ipd(^-d)) is the mean density of 
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impurity for the interacting system. Due to the inter- 
channel coupling, the impurity has a non-zero probability 
TTciosed (studied in \W\ ) to be tightly bound within a closed- 
channel molecule, where it cannot contribute to G(x) and 



to pd- In terms of the probability ir 



1 - ^closed for 



the impurity to be in the open channel, one finds that 
Pd = ftopen/V' I n the thermodynamic limit, the open- 
channel probability is related to the quasiparticle residue 
as 



' d 3 kd 3 q 



(27r) 6 V #q-k + £k ~ £q - A^pol 



(8) 



while the pair correlation function is 



G(x) = 1 



Z 



/0 7T O po 



- 2/(x) 



1 d 3 q 

(27T) 3 

'd 3 k 



(2tt)3 

where we have introduced the functions 



l/k(x)| 



AM = 

/k(x) = 
/(x) = 



d 3 k 



7 A, 



(2tt) 3 £ q _ k + e k - £ q - A£ pol 



' d 3 q 



iqx 



(2tt) 3 £ q _ k + e k - £ q - AE po i 



'd 3 kd 3 



e I ( k -'>l' x /D q 



(2tt) 6 £ q _ k + e k - e q - AE po i 



(9) 

(10) 

(11) 
(12) 



Interestingly, the contribution involving /(x) is an inter- 
ference effect between the subspaces with zero and one 
pair of particle- hole excitations in the Fermi sea. In Fig. [2] 
we plot the numerically obtained deviation G(x) — 1 of the 
pair correlation function from unity (which is both its large 
distance limit and its non-interacting limit): The presence 
of the impurity induces oscillations in the fermionic den- 
sity that are still significant at distances of several l/k F . 

Properties of the pair correlation function. A 

first property of the G(x) function is the sum rule: 



d 3 x[G(x) - 1] = 0, 



(13) 



where the thermodynamic limit was taken and the integral 
is over the whole space. This sum rule follows directly from 
the integral representation of the Dirac delta distribution 
Jd 3 xexp(ik-x) = (27r) 3 5(k). 

A second property is that, in the limit where x — > +oo, 



G(x) - 1 



>+oo 



A4 + B4 cos(2k F x) 
(k F x) 4 



(14) 



The prefactor of l/(k F x) 4 is thus a periodic function of 
x of period Tr/k-p, with a mean value A4 and a cosine 
contribution (of amplitude -B4) reminiscent of the Friedel 
oscillations. The fact that the mean value A4 differs from 
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Fig. 2: (Color online) Deviation of the pair correlation func- 
tion from unity (its uncorrelated value), for M/m = 6.6439, 
= 1 and l/k^a = —2. (a,b) Black thick line: Nu- 
merical solution (J9]). (b) Red solid line: multiscale predic- 
tion (|3ip. Dashed line: non-oscillating bit of the asymptotic 
prediction (|14|l with A4, given by (|26[) . Dotted line: non- 
oscillating bit of the prediction (|30|) . (i), (ii) and (iii): zones of 
the multiscale structure of G(x) defined below Eq. (|3ip . 



zero has an important physical consequence: It shows that 
the polaron is a spatial extended object, since even the first 
moment (x) of G(x) — 1 diverges (logarithmically) in the 
thermodynamic limit. 

Eq. (|14p results from an asymptotic expansion of 
P0I11I12[) in powers of 1/x, obtained by repeated integra- 
tion by parts as in |22j , always integrating the exponential 
function e ik x or e ±lq x to pull out a 1/x factor q 



/k(x). 

/(*)< 



p/{2ir 2 h 2 ) 
x 2 k F T>^ 

-pk F /(2ir 2 h 2 ) 



u=±l 



2 k 2 V 



F(k Fl u;q,u'y kFXU 
J2 F(k,u;k F7 u')e 2kFXU \ 



(15) 
(16) 



=±1 



-p/(8ir 4 h 2 ) ^ 



F(k F ,u;k F ,u r )e lkFx{u ~ u,) .(17) 



Here e z is the unit vector along z, u is the cosine of the an- 
gle between x and k, v! is the cosine of the angle between 
x and q, and the function F is defined as follows: 



F{k,u;q,u') = 

;-M q 2 
, + M~k? 



Am 2 



1 



(m + M) 2 
2p q 



(l-u 2 )(l-u> 2 )^ 



k 2 
2 1 -1/2 



(18) 



Therefore, in Eq. Q the integrals containing |/ q (x)| 2 and 
|/k(x)| 2 provide a 1/x 4 contribution with an oscillating 

2 One also uses the fact that, uniformly in the integration do- 
main: Vrii,n2,n3,ri4 G N, there exists C ni n 2 n 3 n 4 £ K + such that 
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Fig. 3: (Color online) Coefficients A4 and B4 in the asymptotic 
expansion (|14|l of the pair correlation functions, for various 
mass ratios: M/m — 0.1505 (column a), M/m = 1 (column 
b), M/m — 6.6439 (column c), and various Feshbach lengths: 
&F-R* = (row 1), k ¥ R* = 1 (row 2), k F R, = 10 (row 3). 
Upper and lower black lines (solid for Z > 1/2, dashed for 
Z < 1/2): Numerical evaluation of A4 and B4, respectively. 
Red lines: leading order analytical results H26fl and (|27|l . that 
diverge for 7?* > when s of Eq. (|19[1 tends to unity, and for 
R, — when l/(fc?a) — > 0. Vertical dotted lines: polaron-to- 
dimeron crossing point. 



prcfactor, as /(x) does. In Fig. [3] we plot A4 and -B4 
as functions of 1/kpa for various values of the mass ratio 
M/m and reduced Feshbach length k-pR*. 

In the weakly attractive limit. Following 
Ref. [TU], we define for a < 



s = ii (-0^)1/2^ 

m 



(19) 



and we take the limit a — > _ for fixed s < 1 (which 
implies i?* — > +00). Then A£"p i tends to 0~, so that the 
integral appearing in (U) is bounded [10] . we have 



gVq -> l-( S(? /fc F ) J 

a^0~ 



(20) 



and, after integration over q in Eq. ([3]), we get as in [T 

arctanh s 



ft fcf. fc F a 1 

/i 7T s 2 



- 1 



(21) 



Also, the interchannel coupling amplitude A scales as 
(— /cfo) 1//2 , see Eq. ([T]), and a systematic expansion of ob- 
servables may be performed in powers of fcpa, treating 
the interchannel coupling KWud + h.c. within perturba- 
tion theory The terms neglected in the ansatz ([2]) have 
an amplitude 0(A 3 ), that is a probability 0(A 6 ). One can 
thus extract from Eq. © the exact value of Z up to A 4 : 



Z a->0- 



k^a ( kpa \ 
ci— + c 2 — +0(fc F a) 3 

7T \ 7T 



(22) 



The first coefficient of the expansion (|2"2")l is given by 

1 arctanh s 



ci 



1-s 2 



(23) 



it originates from the second term in Eq. (O, a closed- 
channel contribution. This is why c\ = on a broad 
Feshbach resonance where s = 0. The second coefficient is 



C2 
1 

2 



1(1 + a) 2 (s arctanh s — a arctanh a) 
2 1 — s 2 s 2 — a 2 

2m 2 ( arctanh s \1 f 1 + s 2 
+ — ' l] Jl^f- 



arctanh s 



1 s 2 (l — a)(l + a) 2 arctanha 

2(1 -s 2 ) + (l-s 2 ) 2 (s 2 -a 2 ) 

[,s 4 + (1 + a) 2 ,s 2 — a 2 ] arctanh s 
2s(l-s 2 )(s 2 -a 2 ) : 



(24) 



where the mass contrast a = (M — m)/ (M+m) also obeys 
2 arctanha = ln(M /m) 0. The first term in Eq. ((24]) orig- 
inates from the last term in ©, an open channel contri- 
bution. For a broad Feshbach resonance (s = 0), it is non- 
zero, whereas the sum of the other terms of (f2"4"]l . originat- 



ing from the closed channel, vanishes, and C2 — ln ( M / m ) 



It is instructive to analyze the perturbative expansion 
in the exactly solvable limit of M/m — > +00: The im- 
purity can then be considered as a pointlikc scatterer of 
fixed position, for convenience at the center of a spherical 
cavity of arbitrarily large radius R, imposing contact con- 
ditions of scattering length a and effective range — 2i?» on 
the fcrmionic wavefunction (19) . In the thermodynamic 
limit, one can then construct the Fermi sea of exactly 
calculable single-particle eigenstates in this scatterer-plus- 
cavity problem. As shown in [5] on a broad Feshbach res- 
onance, the truncated ansatz ^ provides a good estimate 
of AEp i for M /m — > 00. On the contrary, we find that it 
is qualitatively wrong for the quasiparticle residue: From 
Eq. ((S]), it predicts a non-zero value of whereas the 
exact Z vanishes for M/m 00, which proves the disap- 
pearance of the polaronic character. For an infinite mass 
impurity, indeed, Z is the modulus squared of the overlap 
between the ground state of the free Fermi gas and the 
ground state of the Fermi gas interacting with the scat- 
terer. This overlap was studied in |23j . and vanishes in the 
thermodynamic limit, a phenomenon called the Anderson 
orthogonality catastrophe. Satisfactorily, the perturbative 
expansion ([2"2"]) is able to detect this catastrophe: 



C2 



M/m- 



ln(M/m) 
(1-s 2 ) 2 ' 



(25) 



Such a logarithmic divergence with the mass ratio was al- 
ready encountered in the context of the sudden coupling of 
a Fermi sea to a finite mass impurity, see the unnumbered 
equation below (4.3) in [24]. It originates from the first 



3 Contrarily to a first impression, C2 has a finite limit when s — > a. 
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term of that is from the last term of (O, where it is 
apparent that J d 3 kd 3 q/ (k 2 —q 2 ) 2 diverges logarithmically 
at the Fermi surface. 

Turning back to a finite mass ratio M/m, we now calcu- 
late the coefficients A4 and B4 in Eq. (fH| to leading order 
in kpa. The integrals appearing in Eq. §§§ contribute to 
these coefficients to second order in kpa, see Eqs. (jf 511 61) . 
since 1/2? q scales as grQ. Also Z/iv oven deviates from unity 
only to second order in kpa, see Eq. l[8]). The leading order 
contribution to A4 and B4 is thus given by the interference 
term /(x), and from Eq. (|17|) we obtain 



3s 2 / arctanh; 



B 4 



1 — s 2 \ s 
3 / M\ k F a 



- 1 



2tt 



mil 



2 ' 



, (26) 
(27) 



These analytical predictions are satisfactorily compared to 
the numerical evaluation of A4 and B4 in Fig. [31 

Taking again M /m — > +00, we find a severe divergence 
in the small- k-pa expansion: The coefficient of the term 
linear in kpa in B4 diverges linearly with the mass ratio. 
Instead, the coefficient of the term linear in kpa in A4 
(not given) diverges only logarithmically with the mass 



ratio ' 



t (0) \n(M/m) 

7I-(1-S 2 ) 



This suggests that, for an impurity 
of infinite mass, the asymptotic behaviour of G(x) — 1 
is no longer 0(l/x 4 ) and that it decreases more slowly. 
To confirm this expectation, restricting for simplicity to a 
broad Feshbach resonance (i?„ = 0), we have calculated 
the exact mean fermionic density in presence of a fixed 
pointlike scatterer, obtaining an expression equivalent to 
the one of §2.2.2 of [25] and leading at large distances to 



G(x) - 1 



x-y+oo 2{kpx)'~ 



Re 



2ikpx 



kpa 



1 + ikpa 



(28) 



As this 1/x 3 law has a zero-mean oscillating prcfactor, 
G(x) — 1 has a non-diverging integral over the whole space. 

Multiscale structure of G(x). — In the weakly at- 
tractive limit, it is apparent from Eqs. (|10I11I12[) that the 
functions / q , /k and / are of order one in k-pa, so that the 
leading order contribution to G(x) — 1 in Eq. ([9]) originates 
from the interference term oc /(x) and is also of order one: 



G(x) 



1 



a->0~ 



0(kpa). 



(29) 



One then expects that, when kpx ^> 1, G(x) — 1 drops 
according to the asymptotic 1/ 1 x 4 law with a coefficient of 
order one in kpa. This simple view is however infirmcd 
by Eq. (|26[) . where A4 has a non-zero limit A± for a — > 
0~. This suggests that the 1/x 4 law is only obtained at 
distances that diverge in the weakly attractive limit. 

To confirm this expectation, one calculates /(x) to first 
order in kpa for a fixed x [by using (|2T)j) and neglecting 



4 The integrals over k of the modulus squared of (115ft . and over 
q for the modulus squared of ( 1161 . have contributions scaling as 
(kpa) 2 ln(fcpa) that cancel in their difference. 



AE po \ in the denominator of (|12p]. then one takes the limit 
kpx 3> 1. As shown in the appendix, this gives 



lim 

a->0- 



G(x) - 1 



.4 



(0) 



2M 



7 



30 (kpx) 4 
\n(2kpxM/m 



M 
Am 



cos(2fcp£) 



O 



2(1 - s 2 ) 



(kpx) 5 



(30) 



with the positive quantity e = — AE po \/ Ep <C 1 and 
7 ~ 0.577 215 is Euler's constant. For that order of tak- 
ing limits, the oscillating bit still obeys a 1/x 4 law, with 
the same coefficient B4 as in (|27p : on the contrary, the 
non-oscillating bit obeys a different \nx/x 4 asymptotic 
law (dotted line in Fig. [2p), which shows that the validity 
range of the 1/x 4 law is pushed to infinity when a — > 0~ . 

Remarkably, by keeping AE po \ in the denominator of 
(I12D . one can obtain, see the appendix, an analytical ex- 
pression for G(x) — 1 that contains both the lna;/:E 4 and 
the 1/x 4 laws as limiting cases, and that describes the 
crossover region with cosine- and sine-integral functions: 



G(x) 



A[°\ 



' k F x>l (kpx) 4 



( 77i r 



1 



l ' 2 ) - - Ci(A: F xe) 



f| - Si(fc P xe/2)) sin(fc F a;e/2) 
- Ci{k F xe/2) cos{k F xe/2) +0(1)| (31) 

where the remainder O(l) is a uniformly bounded func- 
tion of kpx > 1 and e <C 1. This formula satisfacto- 
rily reproduces the numerical results, see Fig. [2Jd, where 
e ~ 0.160. It reveals that the pair correlation func- 
tion has a multiscale structure for a weakly attractive 
interaction, with three spatial ranges 0: (i) the logarith- 
mic range, 1 < kpx < e" 1 ^ 2 , (ii) the crossover range, 
(T 1 ! 2 < kpx < 16e _1 , and (iii) the asymptotic range, 
16e _1 < kpx. The logarithmic range is immediately re- 
covered from Ci(u) = lnu + O(l). The e -1 / 2 scaling of 
its upper limit is intuitively recovered if one assumes that 
the relevant wave vectors in (|T2")) obey |k — q| w 1/x: Ne- 
glecting A-Epd with respect to £?k-q in the denominator 
of (fT2]) then indeed requires kpx < (Me/m) 



-1/2 



Conclusion. — The fermionic polaron, composed of 
an impurity particle dressed by the particle-hole excita- 
tions of a Fermi sea close to a broad or narrow Feshbach 
resonance with zero-range interaction, is a spatially ex- 
tended object: The density perturbation induced by the 
impurity in the Fermi gas asymptotically decays as the in- 
verse quartic distance, with a spatially modulated compo- 
nent reminiscent of the Friedel oscillations. In the weakly 
attractive limit, kpa — > 0~ with |a|i?„ fixed, where sys- 
tematic analytical results are obtained, this density per- 
turbation reaches its asymptotic regime over distances di- 
verging as 1 / ( fcp I a I ) and exhibits at intermediate distances 



5 The first relative deviation of (kpx) 4 [G(x) — 1] from its x — > 00 
limit is -24/(fc F xe) 2 + 2^ sin(fc F a:€). For me/M < 6, this is < 10% 
for kpxe > 16. From a similar first-deviation analysis, being in the 
logarithmic range actually requires fcpxe 1 / 2 < (^t/Af) 1 / 2 . 
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a rich multiscale structure. This may have important con- 
sequences on the interaction between polarons |26j . 
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Appendix. — The integral over k and q in /(x) re- 
duces to a triple integral over k, q and the angle 9 between 
k and q. Taking A = k — q| rather than 6 as the vari- 
able, and changing the integration order, we get /(x) = 
(-T/x) / +o ° dAsin(Aa;)^(A) with T = (1 - s 2 )A { 4 0) pe/2, 

pi r^+g 
f(X) = I dq dk- 

J max(l — A,0) J max(A — 



kg /(I - s 2 q 2 ) 
fj-A 2 + k 2 - q 2 + e 



but with a jump J = ip" (2 + ) — ip" (2~) = M l[Am(l-s 2 )]. 

e— >0 

Triple integration by parts over each interval gives 
jr r r+oo 



J cos 2x + ip"(0) + / d\ ip {3) (A) cos Aa; 



/(x) 



The contribution of J reproduces Eq. (j2"T)l . We find that 
<p(A) varies at three scales, e, e 1//2 and e°. For < A < e 3 / 4 , 
we use the scaling A = et and expand (p^ 3 \X) in powers 
of e at fixed t. For e 3 / 4 < A < e 1 / 4 , we use the scaling 
A = e 1 / 2 u and expand ip^ (A) at fixed u. For e 1 / 4 < A < 1, 
we directly expand at fixed A. With 77 = m/M, this gives 



^ (3) (A) = - 



4/(1 -s 2 ) rjY(t-l)/t + 0(2t+l) 



-3 



^ (3) (A) = - 



e 2 {l + 2tf 2e(l- 
(1 - rju 2 )(l + At-ju 2 + n 2 u 4 ) 



s 2 ) 
3 



0(u 4 ) 



^ (3) (A) 



2e 1 / 2 u 3 (l + r 1 u 2 ) 2 (l 



-s 2 ) 4w 4 (l-s 2 ) 
e[l + Q(A)] 
2A 3 (l-s 2 ) 



where Y is the Heaviside function, and the 0( ) apply 
for t G R+, u G M + and A G [0,1], respectively. If one 
needs the integral of (p( 3 \X) cos(Xx) up to an error 0(1) 
uniformly bounded in x and e, one can apply several sim- 
plifications over each interval. E.g., for e 1 / 4 < u < 1, 
one can approximate ip^ 3 \\) by an expansion in pow- 
ers of u, and for 1 < u < e~ 1//4 , by an expansion in 
powers of 1/u. Adding contributions of all intervals, 
we concatenate them by pairs, and further noting that 



J = 0(1), £ 
L 



in 



+°° dt 
-1/2 



0(1), J, 



+ 00 



3,-4 
2 L 



(t + 1/2)- 



0(1), 



G(x) - 1 = 



+ V 



A 



dA|<^ 3 )(A)| = 0(1), we get 
dt cos(xet) 



dX 
2X 



cos(Xx) 



2e (t- 

dX 
2X 



1/2)= 



cos(Aa;) +0(1)}. 



Explicitly evaluatini 
Finally, to obtain 
of (fi(X), so that it is 
{77[ln( 7? A/2)+3/2]-s 
locally split ip" (A) as 
a C°° function. The 
triple integration by 
the singular part, (1 
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